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A mathematical formulation of the overall ship's response 
to bottom impact slamming in regular waves is developed. The 
hydrodynamic problem concerning the definition of the loads 
is first discussed, and a particular physical model is adopted 
for finding the hull vibratory behavior. 

Peasea On the given formulation a general procedure lead-— 
ing to the time history representation of the midship's bending 
moment is suggested. Finally, an illustrative example of 
application to A Mariner ship is described and the results 


compared .with available data. 
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NOMENCLATURE 


ship's sectional area £7 

ship's beam ft 

damping coefficient per unit length ton-sec/ft 
men mode generalized damping ton-sec/ft 
structural damping per unit length ton-sec/ft 


modulus of elasticity ton/£t? 


acceleration of gravity ft/sec? 


mass moment of inertia of hull per unit length with 


respect to an axis normal to the x-y plane; ton-sec* 


sectional area moment of inertia ft? 


hull's shear rigidity ton 


gee mode generalized spring constant ton/ft 


ship's length ft 

length between perpendiculars ft 

added mass per unit length ton-sec*/ft* 
ship's mass per unit length bone gat /Ee 
bending moment ft-ton 

gee mode moment spatial weighting function ton 
total slamming load per unit length ton/ft 
inertial component of P per unit length ton/ft 
buoyancy component of P per ae Wengen econ, fe 


que mode generalized deflection ft 


qc mode generalized forcing function ton 


Mime. vartalple sec 





= | 
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vip Ss UrdaLe £6 

ship's speed ft/sec 

shear force in y direction ton 

ie mode shear spatial weighting function ton/ft 
threshold velocity ft/sec 

regular wave frequency rad/sec 

frequency of encounter rad/sec 

‘qa mode natural frequency rad/sec 

wave length ft 

longitudinal position along the hull ft 

on mode shape, dimensionless 

vertical elastic deflection normal to te 

relative motion between the hull and the wave surface ft 
relative velocity between the hull and the wave surface 
ft/sec 


effective mass per unit length ton-sec* /£t* 


iio mode generalized mass conegac’ (Pe 


component of the slope of y due to bending only rad 
Gensity of sea water lb-sec*/f£t* 


phase angle of z and z motions rad 





fm LNTRODUCTION 

Under certain severe sea conditions the phase relation 
between the bow motion and the surface of the oncoming waves 
is such that an impact may occur between the ship structure 
and the water. This impact of any portion of a moving ship 
upon the surface waves is commonly called "slamming". It is 
in most of the cases a result of large pitch and heave motions 
that force the ship's forward bottom to emerge and re-enter 
the water after hitting its surface. In general this is known 
as"bottom impact slamming’. However, bottom impact is not a 
reguired condition for slamming to occur, and we can talk of 
bow Elare slamming when there is a sudden change of the accel- 
eration of the ship's bow without its actual emergence. 
Finally, as a third type of violent interaction between the 
ship's hull and the seaway we can consider waves breaking 
over the deck, sometimes called "shipping of green water". 

In all cases the impact is such that certain portions 
of the ship's structure generally forward have to sustain 
heavy impulsive loads from the water, possibly causing serious 
local damage along with transient vibrations of the hull. The 
danger of heavy damage or high vibratory response is generally 
reduced at sea by changing course or speed, or both, and this 
means that slamming not only poses a serious threat to safe 
Pier won Olt di SO.resules and considerable Voss i6f time 


that affects the operational capabilities of a ship. 





At the design stage slamming should be studied so that 
its effects on one ship structure and operation can be con- 
veniently estimated. This study could eventually result in 
an improved design, with better handling characteristics, 
better structural reliability and more precise rules for what 
to do at sea in order to find an acceptable.equilibrium between 
safety and effectiveness. 

Essentially we can consider two problems when dealing 
with the ship's response to slamming, a localized one, also 
called a "micro", and an overall or "macro-problem". The 
localized approach deals mainly with local plate forces and 
damage resulting from direct application of the load and it 
tends to be a fairly sophisticated ee aioe problem 
mmen plastic structural analysis playing a possibly important 
role [1].* The overall response involves hull vibrations and 
large midship stresses and bending moment that can be detri- 
mental to the structure as a whole. Following the initial 
localized response and before an overall effect is felt we 
can also define a transition period in which the stress wave 
is travelling along the hull. This logical breakdown of the 
slamming response is useful since it not only makes a study 
approach easier but also tends to establish a clearer boundary 
between the two main sciences involved, hydrodynamics and 


Serictural mechanies. 


*Numbers in brackets [] designate references. 





This study will focus on the overall ship response to 
bottom impact slamming, which is the type more likely to pro- 
duce a real threat to the ship’s structure. The reason for 
treating the problem from the macro point of view is that most 
of the investigations on bottom impact have so far dealt with 
local effects [2], but as ship's dimensions and speed require- 
ments increase the overall response gains a stronger impor- 
mance . 

In our attempt to make a brief review of what we consider 
the most significant papers written on the particular field of 
overall response to bottom impact, we point out first the 
paper by Andrews [3] where a method for determining the elas- 
tic body response of a ship to a seaway was developed. The 
ship was divided into egqual sections and for each one a 
forcing function and mass-elastic parameters representing the 
hull were computed. The forcing function was considered to 
include an unsteady hydrodynamic component obtained from the 
measured rigid body motions, and a hydrostatic component. 

The ship was essentially treated as a free-free beam and the 
waves were assumed to be trochoidal. The damping coefficient 
was given the value zero throughout the computations. 

The work of Leibowitz [4, 5] represents a refined exten- 
Sion of Andrews' report with more precise treatments for 
certain parameters, like added mass, description of the waves 


and damping. The computations were based on an experimental 
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knowledge of the ship's rigid-body motions relative to the 
waves, which were taken from graphical records of periods of 
time where slamming was actually known to have occured. 

The need to represent a ship by a more realistic model 
other than the pure free-free or Timoshenko beam has been 
extensively studied from the vibrations point of view by 
several authors who will be referred to later. Kline and 
Clough [6] specifically attempted the use of an improved 
model to study the dynamic response of the hull, including 
slamming. The influence of bulkhead spacing, machinery and 
Sango location, hull girder and double-bottom stiffness were 
particularly analyzed by these authors. Their model included 
a main elastic axis representing the primary flexural behavior 
of the hull girder and an additional axis accounting for the 
double bottom, which is in fact capable of independent local 
deformations. At intervals the two axes were coupled either 
by rigid connections representing the bulkheads or by springs 
allowing a certain degree of relative movement of the two 
girders. Finally, buoyancy springs provided a support to the 
model. The slamming load or forcing function was assumed to 
be a half sine wave impulse applied to the flexible bottom 
structure at a particular station near the bow. 

The work of Kaplan [7] evaluates existing mathematical 
models describing ship-wave interaction, and it is mostly 


concerned with the computer simulation of wave-induced 





1eiL 


structural loadings, including not only bottom impact slamming 
but also bow flare slamming and springing, or wave-excited 
main hull vibration. The structural model adopted is the 
free-free beam and a definition of the bottom impact forcing 
function similar to the one given by Andrews and Leibowitz is 
cemerowed. The author concludes it 1s not in practice con- 
venient to represent the bending moment due to slamming in a 
spectral form, and suggests as more effective and correct the 
time history representation, which can be directly combined 
with the time history output due to wave-induced effects. The 
computational method requires a hybrid linkage system that 
complements the digital computer. The random wave input is 
obtained through a white noise generator adequately filtered 
to produce a desired spectral form. 

A probabilistic approach to bottom impact slamming has 
been made by Ochi [8, 9, 10] but its analysis does not exactly 
fit the "overall" concept, in the sense used here of hull 
Vibrations and large midship bending moment. Ochi essentially 
finds expressions for the probability density functions of 
occurrence of slamming, of impact pressure applied to the 
ship's forward bottom, and of extreme slamming pressure. These 
results are mainly of interest when studying local bottom 
plate damage, and the overall study as understood here should, 
on the other hand, provide us ultimately with the characteris- 


tics of the midship's bending moment as a random process. 
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Besides the problem of the overall ship response to 
bottom impact slamming, for which we just presented the most 
representative existing investigations, reference should also 
be made to springing or wave-induced main hull vibration, 
Since these two subjects show strong analogies. In fact, 
both seek an overall structural response that involves hull 
wioracions resulting from shi p=waveenoncerceeionr and both 
require a suitable physical modeling of the ship. 

Goodman [11] attempted to find a mathematical method 
capable of solving the springing problem and got results for 
regular waves and irregular seas. The beam model was used 
based on the assumption that the two mode model of vibration 
predominates with a small ee wave energy in 
the higher modes. 

Van Gunsteren [12] develops also a theory for springing 
using a beam model in which the mass and stiffness distribu- 
mioms are discretized. His approach is similar to the one o£ 
Goodman but the numerical results do not completely agree. 

From the considerations of the existing studies on ship 
overall response to bottom slamming we may make some general 
conclusions. First, a general method capable of evaluating the 
response characteristics of a ship in the design stage does 
not effectively exist. Second, the possible theoretical models 
May become extremely complicated if all important parameters 


are to be included and so simplifications are necessary. Third, 
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two main technical areas are involved in the problem, speci- 
fically hydrodynamics and structural mechanics, reflecting 
somehow the complexity already pointed out. 

The purpose of this study is the development of a simple 
and convenient mathematical model capable of effectively 
defining the ship's overall response to bottom impact slamming. 
The desired mathematical model is such that not only new 
designs can be tested as far as slamming performance is con- 
cerned, but also a best design may eventually be achieved. 

The hydrodynamic problem of slamming will first be treated, 
followed by the structural problem of ship hull vibration. 
Finally a numerical example of application will illustrate 


the proposed method. 
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IIT. BOTTOM IMPACT SLAMMING 
meee NECESSARY CONDITIONS 

Several authors have studied and tried to formulate 
mathematically the conditions that are required for bottom 
impact slamming to occur. The leading results we will refer 
to were given by Tick and Ochi. 


Tick {13] established three necessary conditions: 


a) The relative velocity between the bow and the sea 
surface must exceed a critical amount at the time 


Si COntact. 


b) The bow must be out of the water previous to the 


Geneacte 


c) The angle between the wave at some chosen contact 


point and the keel line must be small. 


From tests conducted in regular waves Ochi [14] inferred 


the three following requirements: 
a) bow (forefoot emergence) ; 


b) certain magnitude of vertical relative velocity 


between ship bow and wave; 
c) unfavorable phase between bow and wave motions. 


After conducting tests in irregular waves Ochi could 
reduce the number of his previous necessary conditions by one, 


obtaining the following two: 
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a) bow (forefoot emergence) ; 


b) a certain magnitude of relative velocity between 


wave and ship bow. 


Ochi claims that his first condition is in fact a pre- 
requisite for slamming since tests revealed that slamming 
never occured without bow emergence. This was found to be 
valid no matter what sea state, ship course and speed or 
loading distributions. Being a necessary condition, it is 
not, however, a sufficient one, and a certain magnitude of 
relative velocity between ship bow and wave was found to 
also be required. Ochi called this critical relative velocity 
below which slamming does not occur the threshold velocity, 
which we denote by v* (ft/sec). Tests conducted by this author 
under several conditions with a MARINER model showed that the 
value of the threshold velocity is nearly constant with an 
average of 12 ft/sec for a ship 520 ft long. This value 
should be modified according to the Froude scaling law for 
ships of different lengths. 

Ferdinande [15] further investigated the nature of the 
threshold velocity and found in his tests that it was depen- 
dent on ship speed for different levels of slamming severity. 
The lower the ship speed the higher the threshold velocity 
for the same severity of slamming. This author also Dolieee 
out that this threshold velocity as defined by the Ochi 


probabilistic analysis was not really a strict minimum value 
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of re-entry velocity under which no slamming of the considered 
degree of severity could occur, and he considered it a 

rather artificial concept. Finally he also showed that to 
some extent the squared threshold velocities increased in a 
linear way with increasing slamming severity, related in this 
case to deceleration. 

On his investigations on the circumstances leading to 
slamming, Ferdinande found that a slam could occur after a 
partial emergence of the keel, but with the forefoot being 
already immersed. Also a few slams of a lower degree of 
severity occurred at a moment the whole keel and the forefoot 
seemed to be immersed. The final contact of the bottom with 
Peemvave surface could be located at an appreciable distance 
from the forward perpendicular, but for some slams the exact 
location of impact was difficult to define. 

Another author, Aertssen [16], considers the value of 
the threshold velocity as given by Ochi to be very low, and 
suggests that it should be taken 50%, or 18 ft/sec, greater 
tor the MARINER. 

We may conclude that there is still a certain degree 
of controversy on the exact definition of the conditions 
leading to slamming. Ochi [9] also admits that his threshold 
velocity, which was of course found empirically through exper- 
iments, has not yet been given a rigorous explanation. Anyway 


his two requirements have so far given satisfactory results, 
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and these show a good agreement with practical observations 
related to the probabilistic approach already described. For 
this reason and because they are simple to observe and control 
we will adopt them in our analysis. 

We conclude this section giving a more precise definition 
of the conditions leading to slamming. 

By z(x,t) we denote the relative motion and by 2(x,t) the 
relative velocity between the hull and the wave surface, both 
functions of x, a coordinate defining the position along the 
hull, and t, the time instant. 

The coordinate system used throughout this work = define 
the geometry of the ship is fixed at the calm water level with 
the origin at midships. The x axis coincides with the ship's 
longitudinal axis and its positive direction points to on 
bow. If L is the length of the shap then —b/2) xo & 2. 

The positive z axis points down (this vertical axis will be 
called later y when considering the vibration problem). 

Corresponding to the bow emergence and threshold velocity 
conditions we can write 


ZX) sie) 


Iv 


(la) 


baleen) | v* 


| Vv 


T is the ship's draft at the section where we want to 
investigate if the necessary conditions are or are not satis- 
fied. 

Since bottom impact slamming involves bow emergence, we 


can evaluate (la) at a selected forward location, like the 
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forward perpendicular or station 1, when trying to conclude 
whether the ship is experiencing slams or not. Anyway, atten- 
tion should also be given to other locations along the hull 
towards midships, where conditions (la) can both be met while 
not being so evident at the forward perpendicular. 

Finally we adopt for the threshold velocity the figure 
suggested by Ochi, v* = 12 ft/sec for a ship 520 ft long with 
an Mleopriate Froude scaling law for ships of different dimen- 
Sons, Or 


ve = Omoes (1b) 


2. THE SLAMMING LOAD 

The exact description of the slamming forcing function 
is perhaps the most difficult problem of all this analysis 
and it has not yet been completely explained. Here the dis- 
tinction between localized and overall effects tends to be 
more difficult and a greater degree of simplification and 
assumptions must be followed. Again, most of the investigation 
that has been done has mainly focussed on the local impact 
problem, generally through experimental techniques and empiri- 
cal methods, rather than rigorous ones. 

The first slamming theories [16] were of the wedge impact 
type and were of particular interest for the study of the 
landing of seaplanes. The greatest criticism that can be made 


is that they simply neglected the effects of compressibility 





Ie, 


of both air and water, which tends to become a rather gross 
approximation. 

More recently [17] it has been shown that the air layer 
between the falling body and the surface of the water plays 
a dominant role in the impact phenomenon. The air escape 
from under the body, depending on its shape, is not complete 
and a large bubble provides some kind of cushioning to the 
impact. Under the middle of a flat bottom model the water 
surface is forced to go down by the action of the air, while 
at the same time it rises at the edges, so that the contact 
is actually made before the body cea the undisturbed Sibert 
face. After the impact the situation is more difficult to 
analyze because further air bubbles are produced and more 
realistic three-dimensional tests are hard to control. 

Chuang [18] showed that the cushioning effect of the 
compressible air trapped between the descending body and the 
water surface reduces the maximum impact pressure to about 
one-tenth of the acoustic pressure, and so its presence cannot 
really be ignored when local damage is considered. 

In a more recent paper Chuang '[19] reports the results 
of a series of drop tests with cone-shaped models intended 
to study the three-dimensional effects of slamming. 

The existing studies are somehow difficult to extend to 
the real problem of ship bottom impact, because no account 


1s taken of forward speed or of the three-dimensional flow 
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caused by changes in the shape of the sections. Besides, 

other important factors have been either disregarded or not 
considered in full detail, like shape, velocity and phasing 

of the surface waves, geometry of the hull and hull section 
characteristics, air and water compressibility, two- and three- 
dimensional air and water flow, bubbles and spray. 

We do not know yet for sure which mechanisms are capable 
of playing a dominant role in the impact process and for this 
reason a high degree of simplification must be introduced. 
The factors so far considered tend to be in general important 
from the local point of view and the overall response may be 
simply based on a fairly straightforward physical concept. 
This may be well accepted if the impact process is further 
broken down,as proposed by Oakley [20], into three main phases. 
In phase I the body is approaching the free surface until the 
moment the first contact is made. During this period of time 
the air flow and the wave surface deflection are of predomi- 
nant importance. In phase II the body impacts fully on the 
surface and penetrates it until a more or less wetting is 
achieved. The cushioning effect and the air spray may be 
the prime factors, as well as the water flow around the body. 
Finally in phase III we have the fully wetted problem where 
we can consider the pressures to be static and the forces to 
be the result of the rate of change of momentum. 

HE is difficult to draw the porder between the above 


phases, mainly phases II and III, but we can admit that the 
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ship response is a result of high loads generated as the ship 
re-enters the wave surface and displaces a large amount of 
feet. The physical definition of the forcing function 
becomes rather obvious if we use the concept of added mass. 
and state that part of this force is due to the rate of change 
of the momentum of the added mass of the hull which dives into 
the wave [21]. We will call this the inertial or hydrodynamic 
effect of the slamming load. 

As the ship penetrates the water another force is generated 
because of the increasing buoyancy of the system, and we call 
this the buoyancy or hydrostatic effect component of the slam- 
ming Bead. A more realistic representation of the buoyancy 
term would include the so-called Smith Bee, which accounts 
for the centrifugal forces set up by the orbital motion of 
the particles in the wave, which make the wave pressure not 
conform exactly to the laws of Hydroemment cel 

This model transforms a rather complex problem into a 
fairly simple form which may be incomplete but has the 
advantage of being easily usable in specific numerical appli- 
@eaerons, which is our final aim. tocal effects like air 
cushioning and spray and lateral components of the loads are 
neglected to give a time varying vertical force that has been 
shown to provide acceptable results by Leibowitz [4] and 
Kaplan [7]. 

Reference should again be made to another possible treat- 


ment of this problem as made by Kline and Clough [6], which 
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menas tO be rather oversimplified. In fact, these authors 
simply represent the slamming load by a half sine wave impulse 
applied to the bottom structure at a specified forward station. 
Three durations of loading are considered for the same total 
impulse. 

To conclude the discussion on the bottom impact slamming 
load we will formulate it in mathematical terms so that it 
can be conveniently used to obtain the overall response of 


the ship. 


3. MATHEMATICAL FORMULATION 

As already stated, the total slamming load, which we 
denote by P(x,t) is the sum of an inertial term P, (x,t) and a 
buoyancy or hydrostatic term P, (x,t). These forces are defined 
on a per-unit length basis and they are a function of both 
time and location along the ship's longitudinal axis. 

As defined before, P, (x,t) is the rate of change of the 
momentum of the hull's added mass m(x,t). In mathematical 


Pemmce this can be written as 


P, (x,t) anys (mz) 
3 : 3 ax 
= 2s (mz) + a (mz ) ae (Z) 


imeehne above Gefrini blonetne eee ba aeel  7oueiageuh Dems 
used since the rate of change of momentum occurs not only with 


respect to time but also with respect to x. The added mass 





Ze 


per unit length m(x,t) is calculated for a certain frequency 
of encounter er and it is essentially a function of the shape 
Of the immersed part of the section, which is itself a func- 
Pom of X and t. 

The frequency of encounter Me (rad/sec) is defined in 
terms of the wave frequence w (rad/sec), the ship's speed U 


(ft/sec) and the acceleration of gravity g (Ft/sec*) or 


w. =wt— U (3) 


The term dx/dt in (2) is the relative horizontal velocity 
of the water with respect to the ship. If we neglect the 
particle horizontal velocity, which is in general small, then 
dx/dt % -U, and P, (x,t) becomes 


o 
ox 


3 


are (mz) (4) 


P, (x,t) = (mz) - U 


The buoyancy contribution to the total force per unit 
length, not considering the Smith effect correction, is given 
by 


P. (x,t) = ogA (5) 


© is the density of sea water and A is the immersed cross 
Sectional area of the ship at a particular location x. 

Combining equations (4) and (5) we obtain the final form 
of the total slamming load per unit length for a particular 


location xX: 


g- 
dt 
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Note that in this expression both A and m are functions 
of x and t, since the wave surface position is continuously 


changing relative to the hull. 
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mee, oftP HULL VIBRATION 
ieee HyYSLCAL MODEL 

The problem of ship hull vibration is essentially one 
of finding a physical model that conveniently represents the 
hull's structural behavior and is capable of giving accurate 
or acceptable information about the vibratory response of the 
ship. 

Most of the work on ship hull vibration has treated the 
ship as a free-free beam with the equation governing the motion 
being the Timoshenko equation which incorporates the bending, 
shear and rotary inertia effects. Several numerical methods 
have been used to solve this equation or system of equations 
in order to calculate the natural frequencies and mode shapes, 
and the results agree well with measurements taken aboard 
ships except for the higher modes above the four-node mode, 
where some discrepancies arise. In general, the calculated 
values of the natural frequencies increase faster than the 
measured ones. This-.fact has been explained by several argu- 


ments, mainly in [22]: 


a) inaccurate evaluations of structural parameters 


hike stiffness, load Gdistribucion ocoecammass pete... 


b) vibration has been affected by other factors which 


cannot be accounted for by the simple beam theory; 


c) the assumption that the ship vibrates like a beam 


no longer holds after a certain mode. 
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Hylarides [23] points out shear lag as an important 
factor not taken into consideration by the traditional beam 
approach. 

Also, the beam theory simply ignores the coupling between 
the horizontal and torsional motions and it does not consider 
the influence of local vibrating parts of the hull. Finally, 
we may say that some of the shortcomings of the beam method 
are caused by neglecting the three-dimensional character of the 
ship's structure, like structural discontinuities and three- 
dimensional mass distribution, substituting for it a somehow 
oversimplified physical model. 

Several attempts have been made in order to improve the 
beam modeling of ships for vibrational poeahesd studies. 

Among others, we may refer to the three-dimensional beam-shell- 
Sprung body model of the NS SAVANNAH, with the propeller and 
shafting treated as a sprung body attached to a beam-shell 
system representing the hull [24]. 

Kline and Clough, as already mentioned [6], used two 
beams interconnected by rigid links and springs in order to 
obtain a more realistic model of a ship with a double bottom 
girder capable of independent local deformations. 

The structure of a tanker has also been approximated by 
a model of elastically connected parallel beams [25]. Each 


beam represents the side shells and the longitudinal bulkheads, 
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and the connecting springs represent the transverse bulkheads 


and the transverse rings. 

Hylarides [26] suggests the finite element method as an 
Smecetive means of taking into account the spatial composition 
of the hull, providing therefore a much more realistic three- 
femem>oional representation of a ship structure. With the 
finite element technique, a complex.structure like a ship's 
hull is considered as a composition of elements for which the 
stress-strain relations can easily be formulated. As a result, 
the imperfections of the beam method essentially caused by 
considering the hull as a one-dimensional vibrating beam are 
eliminated. The author also claims that their method also 
takes implicitly into account shear lag effects, local vibra- 
tions and coupling between horizontal and torsional vibrations. 

Based on the finite element technique, a computer program 
called DASH (Dynamic Analysis of Ship Hulls) has been lately 
developed at the Netherlands Ship Model Basin [27]. 

Van Gunsteren [12] compares the beam theory with the 
finite element technique and points out some of the disadvan- 
tages of the latter, mainly from the computational point of 
view. In fact, the finite element method needs three- 
dimensional information about stiffness, wave loads and added 
mass, which is difficult to obtain. For lower modes of vibra- 
tion the finite element technique only gives slightly more 
accurate results but requires considerable computing time 


and this may be its single important disadvantage. 
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From the review of the existing models we conclude that 
the beam theory provides a perfectly acceptable representation 
of a ship structure, at least for the lower modes of vibration, 
which we consider to be the only relevant area for the slamming 
response analysis. The important aspect that influences the 
response is the frequency associated with a particular mode, 
rather than the mode shape itself, since the slamming loads 
have a definite frequency content that is close to the value 
of the first mode natural frequency [28]. 

The beam theory is therefore chosen since it is capable 
of conveniently representing the ship's structure and also 
Mise of the important advantage of providing a straight- 
forward mathematical formulation for the analysis, if certain 


Simplifications can be taken into consideration. 


2. MATHEMATICAL FORMULATION 

The damped vertical response of a ship's hull to tran- 
Sient forces, assuming it behaves like a free-free nonuniform 
beam of length L, eevee by the following system of par- 


tial differential equations [29]: 


yy 
w(x) 2¥ + c(xy + VE) = p(x, t) (7) 
AIS. 
dM (x,t) ay 
sor emma 36170) 92 B36) a (8) 


M(x,t) = EI (x) oy tert) (9) 





EI 
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g t 
oe -— aes + eee) (10) 


distance in the longitudinal direction measured from 
the origin of the coordinate system 

time variable 

vertical elastic deflection, normal to x 

effective mass per unit length, or ship's mass per unit 
length m'(x) plus added mass per unit length m(x). In 
this case we take the added mass as a function of space 
alone based on the ship's calm water waterline sections 
(high-frequency limit) 

damping coefficient per unit length 

shear force in y direction 

total force per unit length due to the ship-wave inter- 
felon, Glven by lo) 

bending moment 

mass moment of inertia of hull per unit length with 
respect to an axis normal to the x-y plane 

component of the slope of y due to bending only 

bending rigidity, where E is the modulus of elasticity 
and I the sectional area moment of inertia 

Shear rigidity, where K is the ratio of the average 
shear stress to the shear stress at the neutral axis 
under vertical load, A is the section area and G the 


shear modulus. 
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The ship is assumed to have free ends so that the boundary 
conditions are: 

Pee /2,t) = V(l72,t) = M(-L72,t) >= M(b/27t)) —e GLaL) 

This system of partial differential equations with its 
boundary conditions has been solved by several numerical 
Meeeods that in general transform them into a set of impliere 
finite-difference equations which can be solved by computer 
[4]. The parameters and variables are then treated as discrete 
rather than distributed. One of the difficulties of the method 
which becomes one of the sources of inaccuracies of the beam 
theory as already pointed out is that it requires an evaluation 
of the yaa dune. and shear rigidity distributions, for which an 
exact procedure is not really available. Besides this fact, 
the solution of a system of partial differential equations, 
which to be of practical use must use a computer, involves 
also approximated numerical methods that are far from accurate. 
However, it happens that the formulation may be considerably 
simplified,as sketched below, in fact leading to a single 
constant coefficients linear second order differential equa- 
tion not involving the parameters EI(x) and KAG(x), if we 
neglect the term involving rotary inertia I. 

The Report of the Committee on Vibration of the Third 
International Ship Structures Congress [22] specifically 
states that investigations showed that there is little effect 


Ge the rotary inertia on the natural frequency of the vertical 
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vibration of destroyers. Van Gunsteren [12] neglects this term 
in his analysis, saying that for prismatic bars its influence 
Smeeme Natural frequency is a quarter of the influence of 
shear, which is already of small influence in comparison with 
bending. 

Accepting then to neglect the influence of rotary inertia, 
the analysis further shows that the dynamic behavior of the 
beam can be presen in terms of series of responses in each 
of its normal modes i, which retain the important property 
of orthogonality with respect to the effective mass per unit 
length [30]: 

ley, 
f U(x) xX. (Ss) XG ee = ee (12) 
ay 2 - J 

Here X, (x) is the normal mode function in arbitrary dimen- 
sionless units, and it simply represents a pattern of relative 
displacements along the length of the beam for a particular 
mode i. 

A generalized coordinate with the dimensions of length 
q, (t) is used to define the displacement time history of the 
system in its ee normal mode. Then the motion in a particu- 
lar mode i is given by multiplying q, (t) by the dimensionless 
normal mode function X (x), and the total response is finally 
given by summing the contributions from all the modes: 


y(x,t) = } 
1= 


qs (t)X, (x) (13) 
al 





By 


Similarly, we can represent M(x,t) and V(x,t) as the 
meaauct of q. (t) by a spatial weighting function M, (x) ona 
Vs (x), respectively, and the form of these functions will be 


determined from the analysis. 


M(x,t) (14) 


lI 
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We assume that P(x,t) can be written in the following 


Semmes form [30]: 


H(x)Q, (t)X, (x) 


E72 (16) 


Piece) =") 


Toe. (Gabe 
1 /D 7 


Multiplying both sides of (16) by X. (x) and integrating 
over the ship's length, the orthogonality property (12) leads 
aman explicit form for the function Q. (t): 

L/2 
OFC ii | P(x,t)X, (x) dx (17) 
-L/2 

Neglecting the term involving IT. and substituting in the 
equations (7) to (10) (where (9) and (10) may be readily com- 
bined) the series representations for y(x,t), M(x,t), V(x,t) 
ame (x,t), (13) to (16), we get the three following equations 
(for simplification of notation we drop here the functional 


dependencies on time and space variables, and we use dots to 


denote differentiations with respect to time): 
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These equations are satisfied if each term in the summa- 


tion is set equal to zero. Combining the resulting equations 





we get 
oe 
2 ad xe 
a : d ~ a dj} V 
pte CAG" 7 Oi SD |e ee Salas 
dx dx 
(21) 
HQ. 

L/2 , 
i UX dx 
-L/2 


If we consider the free motion of the beam with no forcing 
function acting, the right-hand side of (21) becomes zero and 


for a normal mode, after rearranging the equation we get 








ee ods ae a 28s afv)l. 
— 4. °° » &43£°k. 22> ooo ol aerc ies 
qi MAG ax ax 
(22) 
aa 
~_ d°M. 
ne ax 


Since here the left hand side is just a function of time 
and the right hand side just a function of space, we conclude 


that both must be equal to a constant “we? where Ws Saueine 
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natural frequency of the a mode. This leads ae) 








e¢ e 2 — 
SiG am (c/u) Gq. oP ae 0 (23) 
2 2 
aM 2 Gr x 
BY d a d | V | 2 
= EI + EI —/—— =" ae (24) 
dx dx? dx 2 ax RAG) 1. 


Integrating (24) along the length of the beam we get 

x xe 5 
M. = f f Uw X , dxdx | (25) 

-L/2 -L/2 

The following integral gives a more convenient represen- 
Eeeron Lor M.. 

2 
M. = Ww, f (x-s) u(s)X, (s)ds (26) 
-L/2 ; 

pubstituting (24) in (21), multiplying bene sideses 
Xs and taking the space integral of both sides from -L/2 to 
iMemewe get, using again the orthogonality principle (12), 

Becta: eins Pasig Ola vee, 

where Was which we can call the generalized mass, is defined 


by 


i = |i ux, “ax | (28) 
ey : 


and Cas the generalized damping is 


eg. = | ox, “ax 2.9) 
Lo 2n/2 L 


and Ke the generalized spring constant is 


kK. = wot (30) 
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Since the effective mass per unit length yu is only a 

mame tion of x it follows that us 1s a constant for a particular 
mode i and it has dimensions of mass. The generalized damping 
¢, ,which will be discussed in more detail in the next section, 
is for similar reasons also a constant and so the coefficient 
of qs mo (27) is constant. Since Ww, ox the natural frequency 
of the ma mode has a certain fixed value, we conclude that 
(27) is a simple constant coefficients linear second order 
differential equation where the unknown is q, (t) and the forc= 
mieecunction is Q. (t). The solution of such an eguation 1s 
simple if initial conditions are given. At this point we can 


assume that at time t=0 the beam is at rest so that q, (0) = 





q; (0) = 0. The solution is then given in a closed form as 
t Q. (tT) = | 
q, (t) = ! = exp [= (¢,/2u;) (t-t) [sin[A, (t-t) fat (3a) 
ie 
| 2 Pe = 2 
\, =V WwW; (1/4) (c, /u;) (32) 


Knowing q, (t) as well as the normal mode shapes and 
natural frequencies, it 1s possible to compute M. (x) from (26) 
and finally obtain the bending moment from (14) for any loca- 
tion x along the ship. 

The integral definition of M, (x) Given an (26) can be 
interpreted as if we considered the main contribution to the 
bending moment to be a result of the inertial loads distribution 
miong the hull due to the vibratory deflection, ancluding the 


fluid inertial force associated with the added mass. 
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The differential equation (27) and its closed form solu- 
tion (31) represent, then, the analytical solution to the ship 
structural model and now the great advantage of its final 
computational simplicity becomes more evident. The analysis 
will be complete after giving a more rigorous definition to 
the several parameters involved, particularly the added mass 


and the damping coefficient. 


3. DEFINITION OF PARAMETERS 

The solution for the transient response of the ship 
structure to a bottom impact slamming load, as given by the 
normal approach to the beam theory, requires the knowledge of 
a certain number of parameters, namely ship mass and added 
mass distribution, damping coefficient and normal mode shapes 
and natural frequencies. 

The ship's mass or weight distribution may be determined 
through classical methods from a knowledge of its structural 
Main characteristics and internal arrangements. In general, 
this calculation is sae of the design process and we may 
consider it here as a readily available piece of information. 

The added mass distribution is a more involved problem 
and the approach followed here is heavily based on the one 
adopted by the Computer Seakeeping Program developed at M.1.T. 
[31]. Here the computation of the sectional added mass is 
based on the hydrodynamic problem of determining it for an 


infinite length cylinder floating at an infinite depth of 





oF 


water and oscillating vertically with a small harmonic motion, 
the viscous and end effects being neglected. The sections 

are defined as Lewis forms with main parameters area, draft 
and half beam, and the numerical procedure for the actual 
computation of the added mass is the one due to Grim. The two 
main restrictions of this algorithm are taken into account 

in the sense that it cannot be successfully used either for 
very small sectional areas or for bulbous sections with very 
large sectional areas. The added mass resulting from this 
method is frequency dependent (frequency of encounter), which 
should be taken into account when SoinsnelnG the forcing Pine 
miromeP (x,t) as given by (6). 

When seeking the coefficient of q. (t) in (27) the high 
Frequency asymptotic limit of the added mass is used, since 
the ship's natural frequency is relatively high and therefore 
we can assume the added mass to have already reached its 
frequency independent asymptotic value. 

In the next section describing the proposed general pro- 
cedure it will be mentioned how the calculations will be 
Sierred out for specific numerical applications. 

The problem of defining the damping associated with hull 
vibration is quite difficult since its mechanisms have not 
yet been completely explained. 

Im {22] ship hull damping is classified as a combination 


O£ the energy dissipation effects that occur not only in the 
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fluid medium but also within its own structure. The internal 

or structural damping is the result of such varied mechanisms 

as working and fraying of overlapping plates, motion of 

loose cargo in the hull and damping between rivets and struc- 

tural joints. On the other hand, the hydrodynamic or external 
damping can result from energy losses with low frequency sur- 

face wave generation, water friction or losses in the boundary 
layer and high frequency acoustic energy dissipation. 

Investigations have concluded that the hydrodynamic damp- 
ing is less important than the internal damping by at least 
one order of magnitude but the contribution due to the genera- 
e10On MP fue face waves in general may be significant. Also, the 
internal material damping is negligible and cargo damping has 
not yet been really described. Experimental measurements 
give results for the damping coefficient that can be assumed 
to include the effects of the structure, cargo and minor hydro- 
dynamic effects like water friction and the generation of 
pressure waves. 

Goodman [11] uses data from Johnson and Ayling [32] 
Giving in a graphical form the damping ratio due to effects 
other than the generation of surface waves and forward speed, 
as a function of natural frequency. In a discussion of the 
same paper [11], Belgova criticizes this approach, stating 
that the dependence of damping upon frequency should only be 


established from data for a particular mode at a time, since 
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ge CONGitions for the behavior of the ship structure in 
aqifferent vibration modes are not the same. 

In spite of its shortcomings the approach suggested by 
Goodman for the damping coefficient definition will be adopted, 
Since it is to our knowledge the most recent work in which 
this problem is raised. 

We will consider three main factors in the damping coeffi- 
cient: first, a structural damping factor c./u which includes 
also the minor hydrodynamic effects of water friction and 
generation of pressure waves, and which is only a function of 
the mode frequency; second, the generation of surface waves 
effect; and finally, a forward speed correction factor that 
takes the form -U(dm/dx). Note that Co and -U(dm/dx) are 
given on a per-length basis. 

At the relatively high ship's structural mode frequencies 
the effect of damping due to wave generation may be considered 
negligible as compared with the other mechanisms that are 
iaeluded in co/U. In fact, as already mentioned, the struc- 
tural damping may be about an order of magnitude more important 
than the hydrodynamic damping. Then we are left with the 
ve contributions co /u and -U(dm/dx) as defined before, and 


using the definitions (28) and (29) we get 


Cc 
seo, -0 Me ely a (och) 
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_ Heya S 
Cc. = J Ie), - U any x. “ax (34) 





-L/2 
C 72 
oe So age dm a 
ee ia ir, uf ae eee ox (35) 


To complete the knowledge of the parameters required to 
solve completely the vibratory response of the ship, we finally 
need the normal mode shapes X. (x) and natural frequencies W.- 

Several procedures have been developed in order to compute 
X. (x) and War like the Prohl-Myklestad digital method reported 
meme | Or the finite element program DASH already mentioned 
evi. 

The literature reports mode vibration data for various 
types of ships under different load conditions, and so we can 
use it for similar designs if the computer methods are not 
readily available. 

Summarizing this section we may say that the required 
parameters for the ship's transient response solution are 


obtained in the following way: 


a) the ship mass distribution by the standard naval 


architecture procedure; 


b) the added mass using a computer program applying 
the Lewis form concept in conjunction with the 


computation method of Grim; 


c) the structural damping coefficient c./U will be 


Gowen srom a Ghanhical) cepresentaclon -qivingsres 





d) 


4l 


relation to the natural frequency [11]; the general- 
ized damping will be computed in accordance with 


51); 


the mode shape functions and natural frequencies will 
be either computed using available computer pro- 


grams Or derived from existing data for simi lanesiiieee 
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fy GENERAL PROCEDURE 
1. INTRODUCTION 

Having discussed the theoretical aspects of the proposed 
method for the evaluation of a ship response to bottom impact 
slamming, we must now formulate it in a general procedure 
able to be conveniently applied to specific cases. This will 
mainly be a discussion of the sequence of operations to be 
performed, since the mathematical formulation has already been 
given. The method will be further illustrated by a numerical 
application described in Part V. 

Before proceeding, a point to take into consideration is 
the description of the sea conditions for which the response 
is to be evaluated. The most realistic situation and the one 
more likely to give more useful information is the one of 
irregular seas characterized by a certain sea state or energy 
Spectrum. However, the first step to be taken must concern 
the rather simplified situation of regular waves, from which 
a more general irregular-waves result can eventually be 
derived. 

The problem of extending results from regular to irregular 
waves may be quite involved and several approaches may have to 
be tried. In any case, the possibility of a linear dependence 
on the amplitude of the waves should be studied, since if 
we admitted linearity to exist, the use of the superposition 


principle would be possible. In this case, we would consider 





43 


the seaway to be composed of many regular wave components, 

each having its own amplitude, length and direction of travel. 
Response values per foot of amplitude for several frequencies 
in conjunction with a sea energy spectrum would then lead the 
corresponding response spectrum. The acceptance of this 
approach has, however, to be verified in terms of the linearity 
of the response with the wave amplitude, which may not be 
obvious. 

Another possible way, at least formally, would be to 
study the feasibility of treating the ship as a system of 
one degree of freedom acted by a driving force Q,(t), as sug- 
gested by equation (27). Then using spectral analysis con- 
cepts the knowledge of the frequency spectrum of the forcing 
function would lead to the frequency spectrum of the response, 
in this case q, (t). The problem would then be how to find 
Ene frequency spectrum of q, (t) or P(x,t), and hew tower 
from the spectrum of q, (t) the frequency spectrum of the bend- 
ing moment, which is ultimately what we are interested in. 

The present work will only deal with the response in 
regular waves of a given frequency and amplitude, and this 
will hopefully become the basis for the generalization of 
the method in order to handle the more realistic case of 
irregular waves. 

An important point of the Sroeeanude. Ae the computation 


of the ship motions in regular waves, which will be made using 
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the Seakeeping Program already mentioned [31]. Here the pre- 
diction of seakeeping performance is still limited to two 
modes of ship motions (pitching and heaving) for head long- 
Mecteca seas. The theory adopted in the report is essentially 
the linear strip theory of Korvin-Kroukovsky where surge, 
Sway, yaw and roll have been neglected. 

The consideration of only two degrees of freedom in head 
Seas 1s acceptable for the study of slamming since these con- 


ditions lead to the most severe bottom impacts [8]. 


PD EoCREPTION OF STEPS 
The steps involved in the computational procedure are 
given below. For ease of reference they are denoted by the 


capital letters (A) to (G). 


eee nip description 


The following are assumed to be given: main dimensions 
and characteristics, mass distribution, natural mode shapes 
and frequencies, table of offsets or ship's lines and Bonjean 
curves {sectional area curves). 

The structural damping coefficient c./u is taken from 
Figure 2 in [11] for the different natural frequencies. 

The threshold velocity v* comes from (ls) and the nececs 


Sary conditions for slamming are found from (la). 





45 


(B) Regular waves description (frequency) 


The frequency w is selected within the frequency range 
of a sea energy spectrum of particular interest. 


ine frequency of encounter We LS Gerermned sirens 


(C) Added mass 

The frequency-dependent added mass is calculated using a 
computer program we call ADMASS which is based on [31]. 

For each station the computations are repeated for 
several drafts and for each section the required inputs are 
the frequency of encounter, the draft, and the half breadth 
and sectional area known from the ship's lines and Bonjean 
curves. Then for each station we shall have for a certain 
frequency of SMeoWMee 2 set of values of the added mass, each 
corresponding to a particular draft, so that a curve may be 
drawn representing the variation of m withT (Figure l). 

The high frequency limit of m is evaluated using the 
same computer program ADMASS, for every station and for the 
inputs half eer sa draft and sectional area at the calm 
water waterline. The frequency input is increased so that a 
limiting asymptotic value may be obtained. This high fre- 
quency limit added mass will then be taken as being only a 
function of space or m(x), and it will be summed to the ship's 
mass distribution m'(x) in order to get the effective mass 


u(x). The procedure is sketched in Figure 2. 
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Figure 1 Added mass as a function of the draft for station n 
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Figure 2 Effective mass computation 
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(D) Ship motions in regular waves 


The Seakeeping Program [31] gives for each ship station 
the amplitudes and phase angles of the relative motion z(x,t) 
and the relative velocity z(x,t) between the hull and the 
wave. 

mne frequency of the relative motions ange vellocieic ee. 
the frequency of encounter known from step (B). 

The wave frequency also known from (B) 1s a required 
input for the program, as well as a set of offsets which may 
be taken from the ship's lines. If the offsets are not avail- 
able the program may use simply the principal ship's charac- 
teristics. 

For a unit wave amplitude the program gives non-dimen- 
Sional motion response amplitudes, so that an amplitude may 
be selected such that the conditions for slamming (1) are 
satisfied. This means we know the relative motions and velo- 
Cities (amplitude and phase angle) at every ship station when 


slamming 1s occurring. 


foe Loading function P (x,t) 


Knowing the relative motion and relative velocity at each 
station in terms of frequency, amplitude and phase, and since 
these are sinusoidal functions we can generate their time 
mistory. This méans we can have for different equally spaced 
imetants Of time, which we must select, aevalWe lob eho eclae 


tive motion and the relative velocity at each station. 
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mie relative Motion In Conjunction withwthne drabeawil 
give the instantaneous draft T(t), or the immersed portion 
of the hull at a certain station and for a particular time 
instant. 

From (C) we can take the value of the added mass m for 
each instantaneous draft since we have a graphical represen- 
meron OL the variation of m with T. Similarly, from the 
Bonjean curves we can take the value of the sectional area 
For each station and for each instantaneous draft. 

We finally have for each station and for equally spaced 
instants of time the values of the relative velocity 2(x,t), 
the added mass m(x,t) and the sectional area A(x,t). After 
performing the product of the added mass by the relative 
WVetocity (mz) for each time and station, numerical differen- 
tiations are required to finally obtain P(x,t) as given by 
equation (6). 

The procedure is sketched in Figure 3, where I denotes 
the operation of finding the instantaneous immersion and II 


mame Solution of (6). 


(F) Q,(t), M,; (0) and generalized coefficients 
The petsice atin On Q.(t), M,(o), We and k, follow imme- 
diately from equations (17), (26), (28) and (30) respectively, 
the first three involving an integration. | 
The computation of Cas which involves an integration and 


a differentiation, is now possible since we know all the t 


ee 35), 








Figure 3 
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station n 


Computation of the loading function for station n 
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(G) Closed solution of the governing equation (31) 
' and computation of the midships bending moment (14) 


Equations (31) and (14) can be solved in one step, leading to 
ere final result of interest M(0,t), or the midships bending 


moment time history for a particular type of wave. 
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See NUMERICAL EXAMPLE: — APPLICATION TO A MARINER SHIP 

The main results of a numerical application to the MARINER 
femelebe given in this section. 

We should notice it is not worthwhile to publish in 
detail all the calculations, since most of them are of no 
particular significance. Then, besides reporting the more 
important results, the decisions that were taken along the 
procedure will also be discussed and justified. 

The MARINER type was chosen since there are many of these 
ships in operation today, and also in has been extensively 
studied in the past. Examples are the already-mentioned works 


by Ochi [8] and Lewison[17]. 


1. GENERAL PROCEDURE AND MAIN RESULTS 
fay Ship description 

The MARINER principal characteristics are listed in 
Table l. 

The selected speed was 16 knots. For the MARINER the 
speed range of 15 to 20 knots has been considered as critical 
for slamming [33], and this is the reason why the value of 
16 knots was chosen. 

The weight distribution is given in Figure 4. 

The first mode natural frequency is given by Kaplan [7] 
for a MARINER, as CEs he 9.42 rad/sec. This author takes one 


first mode shape to be a simple parabolic curve. 





Table I 


General Characteristics of the MARINER 


Length overall, ft-in 

Length between perpendiculars, ft-in 
Beam, maximum, molded, ft-in 

Depth, main deck at side, ft-in 
Draft, ft-in 

Displacement, ton 

Pekeek coefficient 

mersmatic coefficient 

Midshipsection coefficient 

Number of stations 


Station spacing, ft-in 


oe 
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Figure 4 Mariner weight distribution 
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Data for the vibration modes 2 and 3 was not readily 
available, so that we decided to use the results found for 
a similar ship. Since these modes have a rather secondary 
contribution to the response, at least when compared with mode 
1, the imprecisions introduced should not be very significant. 

The report on the NS SAVANNAH [24] became the source of 
the required information about mode shapes and natural fre- 
quencies. In fact, the NS SAVANNAH and the MARINER types 
are similar in their principal dimensions, as shown in [34]. 
Besides, the first mode natural frequency for the NS SAVANNAH 
1s given as 9.498 rad/sec, which is very close to the figure 
of 9.42 rad/sec mentioned above. For these reasons we decided 
to adopt the data given in [24], as shown in Figure 5. 

The structural damping coefficient c./U for each mode 
natural frequency iS given in Table II. 

The threshold velocity v* is 12 ft/sec, so that the con- 
ditions for slamming may be written as |z| i PI 128 eha\cl 


[2 | fee Lt / Sec. 


(B) Regular waves description 


The selected frequencies and correspondent frequencies 


of encounter are given in Table III. 


(C) Added mass 
The added inass per unit length plotted as a function of 


the draft for the particular case of w = 0.8 rad/sec 





56 


24.609 rad/sec 


9.4968 mad7cee 
= 17.278 rad/sec 
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Structural Damping Coefficient 


mode We (rad/sec) c./u ( /sec) 
IE 9.498 0205s 
2 Vera rs: 0.130 
3 24.609 0.260 
Table IIIf 


Regular Waves Frequencies 


w (rad/sec) oe (rad/sec) 
OZ 0272556 
0.4 0.5344 
O56 0.9024 


0.8 mee s 
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(we aml. s3/5 rad/sec) is given in Pigure 6 ferea selected 
number of stations. 
The high frequency limit of the added mass, as used to 


Sompute the effective mass, is plotted in Figure 7. 


(D) Ship motions in regular waves 


A wave amplitude of 15 ft was found to satisfy the slam- 
Ming requirements defined in (A). To illustrate this fact we 
give in Table IV the relative motion for the first seven 
stations, for the wave frequency of 0.4 rad/sec (wave ampli- 
tude 15 ft). We can see that the requirements |z| 22g ee 
and [Z| > 12 ft/sec are both met for Stations 1 through 4. 
For the remaining stations the amplitudes of relative motion 
and velocity decrease, until they reach a minimum at Station 


Oe. 


fee boading function P (x,t) 

First a simple computer program generates the values of 
Pe mimmnersions at each of the 20 stations, for time intervals 
of 0.2 seconds starting at t = 0.0 seconds. The input is 
the frequency We plus the relative motion amplitudes and 
phase angles as obtained from (D). 

The computer essentially prints out the values of a 
Sinusoidal function of time to which it adds the ship's draft 

ne sin(w.t + G) + 27 


Lor t,=0.0, t,=0.2, t,=0.4, etc. 





HN 
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Figure 7 High frequency limit of the added mass 
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Table IV 


Sample of Relative Motions and Velocities 


V = IG knots (27.0222 fe/cee) 
WL/LBP = 2.4297 

w = 0.4000 rad/sec 

ae 0.5344 rad/sec 

wave amplitude = 15 ft 


Relative Motion Re late veny cloc may, 
Siac i1on Amplitude (ft) Api tude mee, see) 
Ji | 39.08 20.88 
2 Bo. 29 18.86 
3 Sub al 16.89 
4 28.04 14.98 
5 24.61<27 aes 
6 pee S427 11.40<12 
7 18.24<27 See sk 
Table V 
Ei =i, 4 sre 
Mode Prone eCes ee) (tone eee aa (TonvyLe) (ton) 
1 362.45 21.377 32.7 x 10> «1,209 2 To 
2 403.10 54.909 ioe o> ero oe 
3 3 


3 754.09 200.333 450 x 10 -1,568 x 10 
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This first printout gives the instantaneous immersions 
and for each of these we must tabulate the added mass and the 
sectional area. 

A second computer program then performs the numerical 
differentiations and other operations required to obtain 
P(x,t) as defined in equation (6). The differentiation sub- 
routine uses a Lagrangian interpolation polynomial of degree 
4 relevant to five equidistantly-spaced argument values. The 
input includes the frequency Wee the relative velocity ampli- 
maces and phase angles for the 20 stations, and the values 
of added mass and sectional area tabulated as explained before. 

These two computer programs are called I and II respec- 


tively, in Figure 3, where the procedure is sketched. 


(F) Q.{t) | M, (0) and generalized coefficients 

The same computer program II described before uses an 
integration subroutine to compute Q, (t). In addition to the 
defined input it must also read the mode shapes X, (x). 

A third computer program which we call III is just a 
combination of numerical differentiation and integration sub- 
routines, and it gives the values of Ma Cas K. and M, (0) 
listed in Table V. 

The integration algorithm follows Sympson's Rule together 


with Newton's 3/8 rule. 
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(G) Closed solution of the governing equation (31) and 
computation of the midships bending moment (14) 


A computer program which we call IV solves equations (31) 
and (14). The numerical integration method is the one already 
described, and the input includes the results of step (F) 
along with the mode shapes X. (X). 

Obviously II, III and IV may be combined in a single 
computer program as will be referred later. 

The final results we were seeking, or the midship 
bending moment time history for the four input regular wave 


Peequencies, are given in Figures 8 to ll. 


Zee LSCUSSION 

The midship bending moment response to slamming we obtained 
is a time-varying function that presents a periodicity similar 
to the one of the waves, as would be expected. This is parti- 
cularly evident for wave frequencies of 0.6 rad/sec and 0.8 
rad/sec, where the plotted time histories are given for well 
over a period (Figures 10 and 11). 

Referring to Figures 8 through ll we may conclude the 
following: 


a) For w=0.4 rad/sec the bending moment varies between 
-—/Q x 10° ft-ton and zero. 
b) For w=0.4 rad/sec it varies between approximately 


WHEO 4 M0? “EELeGn and COM ORE Ce eone 
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Figure 8 Midship bending moment for regular waves with 
w = 0.2 rad/sec 
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Figure 9 Midship bending moment for regular waves with 
w = 0.4 rad/sec. The contributions of modes 2 and 
3 are also shown. 
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Figure 10 Midship bending moment for regular waves with 
w = 0.6 rad/sec 
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Figure 11 Midship bending moment for regular waves with 
w = 0.8 rad/sec 
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c) For w=0.6 rad/sec the variation is between approxi- 
mately -500 x 10° £E-COM sane ee ua Pe. om. 

fmerinally, for w=0.8 rad/sec the minimom value ic 
-160 x 10° ft-ton and the maximum is 130 x 107 


ft-ton. 


These figures show a strong correlation between bending 
moment amplitude and frequency. The amplitude is a minimum 
for w=0.2 rad/sec, it is a maximum for w=0.6 rad/sec and it 
is of about the same order of magnitude for the remaining 
two frequencies. 

Towing tank tests also show this kind of correlation for 
the MARINER [8, 33]. The most severe slamming was found to 
occur for values of wave length (WL) over ship length (L) 
of 1.0 to 1.25, and almost all the waves of WL/L = 0.75 to 
2.25 and of heights greater than 20 ft were found to cause 
slamming. 


In terms of WL/L we have for our four frequencies the 


following: 
w (xrad/sec) WL/L 
OZ 2p, 7 LENG 
0.4 Sea 297 
Orc IL OY Sis, 


Or 0.6074 
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It is interesting to note that our results agree with 
the towing tank experiments reported above at least in two 


merits: 


a) The response is maximum for WL/L = 1.0799 (w = 0.6 


mad/sec) which lies in the “vange vor (ete eome 2 oe 


b) The response is minimum for WL/L = 9.7190 (w = 0.2 


rad/sec) which lies outside the range 0.75 to 2.25. 


Concerning the relative contributions of the different 
modes to the response we can conclude the following (as an 
illustration the three modes are shown for the particular 


case of w = 0.4 rad/sec in Figure 9): 


a) For w = 0.2 rad/sec the second and third mode 
amplitudes can reach about 20% of the first mode 


amplitude. 


b) For the remaining studiec frequencies the importance 
of modes 2 and 3 decreases but their amplitudes 
still reach a maximum of about 10% of mode one's 


amplitude. 


c) Neglecting the influence of modes 2 and 3 is prema- 
ture, since their contribution bostemcmer aa 


response can still be significant. 


Comparing our results with the experimental ship bending 


mements at midships for a 368 £t donq dectro ee watea ceceume: 
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17 knots [4], we can conclude they are of the same order of 
Heratude. In fact [4] gives a time history similar in shape 
to what we obtained, and the maximum value it reaches is 


3 ft-tons. We should note that this 


approximately 50 x 10 
experimental response refers to a real situation of irregular 
seas (no Sea State given), where the ship's rigid body motions 
were taken from actual records. 

Kaplan [7] obtains through simulation a time varying 
bending moment for the MARINER in irregular seas that seems 
to decay exponentially with time rather than having a certain 
periodicity. The maximum amplitude it reaches is about 


3 ft-tons (for a speed of 12 knots vandesea cea cem yn 


25 Same 
which seems to be guite low compared with our values. 

It is also possible to compare our results with the stan- 
dard bending moment calculation. For the MARINER [35] gives 
these for a wave height of 1.1/L, which is close to the figure 
of 15 ft amplitude used here. The most critical condition is 
for hogging where the value of the midship bending moment for 
faewtnlt-load condition is 388,300 ft-ton. We see then thet 
the maximum midship bending moment due to slamming in regular 
waves may be larger than the maximum standard bending moment 
by as much as about 30% or more, for a particular wave fre- 
quency (in our case for w = 0.6 rad/sec). 


We may conclude that the numerical results of application 


to the MARINER as described in this section are within an 





Wal 


acceptable order of magnitude. However, a more extensive 
comparison with experimental data should be carried out if 


more accurate conclusions were to be taken. 
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VI. CONCLUDING REMARKS 

The numerical application to a MARINER of the midship 
slamming response method outlined in Sections I to IV showed, 
as already discussed, acceptable results within the ranges 
reported in the literature. It is necessary, however, to 
examine again the procedure and study the possibilities of 
expanding and refining it. An important point is weighting 
the validity of the several implied assumptions which include 


erm iy : 


ome the definition of the Meadimg funeeiemeC we cs 


given in (6) with the Smith effect neglected; 


b) the neglect of localized hydrodynamic effects, as 


alr cushioning and spray; 


c) the consideration of only vertical plane motions 
in head seas and the adoption of the Korvin- 


Kroukovsky equations; 


d) the use of the beam theory normal mode of approach, 
with neglect of such problems as local effects, 
coupling of vibration modes and higher modes of 


vibration. 


A desirable extension of the analysis should provide 
for the solution of the irregular seas condition as already 
Seated, since this is the onewilikel Se qavemicone Username 


GBeamustic results. 
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The present work with all its imperfections arrived at a 
formulation and procedure that have the virtue of not being 
difficult to apply to specific cases, and this must be con- 
sidered when estimating the weight of the above assumptions. 

The computational sequence is in general easy to carry 
except for the preparation of data (m and A as a function of 
time for each station) for the so-called computer program II. 
Anybody with a fairly good programming experience may easily 
devise a quicker way of doing it. In fact programs II and 
III are extremely simple and it is possible to combine them 
with the ADMASS program avoiding then the lengthy intérmediate 
tabulation and graphical representation of m as done in 
Figure 6. Also the tabulation of A (sectional area) may be 
avoided by some convenient description of the hull form, as 
done for example in subroutine OFFSET of [31]. The reason 
this was not carried out here is that our main objective was 
the discussion of the method itself and its formulation, 
rather than the preparation of a computational algorithm. 

To conclude, we sincerely hope that in the near future 
the study of this subject is further developed and improved 
Bortchat eventually the response Of Ships) to any type of Slam- 
ming may be readily obtained. If so, the benefits to be 
derived from this analysis would be enormous, with ship's 
structural design certainly closer to an optimum and safety 


and effectiveness undoubtedly increased. 
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